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Abstract: 

In statistical QCD, the thermal properties of the quark-gluon plasma can be determined 
by studying the in-medium behaviour of heavy quark bound states. The results can be 
applied to quarkonium production in high energy nuclear collisions, if these indeed form 
a fully equilibrated QGP. Modifications could arise if an initial charm excess persists in 
the collision evolution and causes quarkonium regeneration at hadronization. 
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1. Introduction 

We know from statistical QCD that strongly interacting matter undergoes a deconfine- 
ment transition to a new state, the quark-gluon plasma. How can we study this state - 
which phenomena provide us with information about its thermal properties? The main 
probes considered so far are 

• e-m signals (real or virtual photons) 

• heavy flavours and quarkonia (QQ pairs) 

• jets (energetic partons) 

The ultimate aim must be to carry out ab initio calculations of the in-medium behaviour 
of these probes in finite temperature QCD. 

In high energy nuclear collisions, we want to study the deconfinement transition and the 
QGP in the laboratory. The ultimate aim here must be to show that experimental results 
confirm the predictions of statistical QCD, or that they disagree with them. If the latter 
should happen, we can unfortunately not conclude that statistical QCD is wrong; a more 
likely conclusion would be that nuclear collisions do not produce the medium studied in 
equilibrium QCD thermodynamics. 

I want to consider here a specific case study for the program just outlined: the spectral 
analysis of quarkonia in a hot QGP and its application to nuclear collisions. 

The theoretical basis for this analysis is: 

• The QGP consists of deconfined colour charges, so that the binding of a QQ pair is 
subject to the effects of colour screening. 

• The screening radius ru(T) decreases with temperature T. 

• When ro(T) falls below the binding radius of a QQ state i, the Q and the Q can 
no longer bind, so that quarkonium i cannot exist pQ. 

• The quarkonium dissociation points T i; specified through ro(Tj) ~ r*, thus deter- 
mine the temperature of the QGP, as schematically illustrated in Fig. ^ 
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Figure 1: Quarkonium spectral lines as thermometer 

Experimentally, quarkonium studies also provide a great tool: 

• In AA collisions, quarkonium production can be measured as function of collision 
energy, centrality, transverse momentum, and A. 

• The onset of (anomalous) suppression for the different quarkonium states can be 
determined and correlated to thermodynamic variables, such as the temperature or 
the energy density. 
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Figure 2: Charmonium survival probabilities vs. energy density 

• The resulting thresholds in the survival probabilities <S$ of states i can then be 
compared to the relevant QCD predictions, as illustrated in Fig. El 

In this way we can, at least in principle, obtain a direct comparison between experimental 
results and quantitative predictions from finite temperature QCD. 

2. In -Medium Behaviour of Quarkonia: Theory 

We consider as quarkonia those bound states of heavy quarks which are stable under 
strong decay; they are thus pairs of charm (m c ~ 1.3 GeV) or beauty (ra& — 4.7 GeV) 
quarks, whose overall masses fall below the open charm or beauty thresholds. The large 
quark mass allows spectroscopy based on non-relativistic potential theory |2j. Hence the 
Schrodinger equation 

Um c - — V 2 + V» W(r) = M&iir), (1) 
I m c J 

using the "Cornell" form for the confining potential j3] , 

V(r) = ar--, (2) 
r 

in terms of the string tension o ~ 0.2 GeV 2 and the gauge coupling a ~ vr/12, determines 
the masses M, and the radii of the different charmonium and bottomonium states. The 
results are summarized in Table 1 and are seen to give a good account of quarkonium 
spectroscopy, with an error of less than 1% in the mass determination AM for all (spin- 
averaged) states. 



state 


JH> 


Xc 


v/ 


T 


Xb 




Xb 


'VV/ 


mass [GeV] 


3.10 


3.53 


3.68 


9.46 


9.99 


10.02 


10.26 


10.36 


AE [GeV] 


0.64 


0.20 


0.05 


1.10 


0.67 


0.54 


0.31 


0.20 


AM [GeV] 


0.02 


-0.03 


0.03 


0.06 


-0.06 


-0.06 


-0.08 


-0.07 


radius [fm] 


0.25 


0.36 


0.45 


0.14 


0.22 


0.28 


0.34 


0.39 



Table 1: Quarkonium spectroscopy in non-relativistic potential theory 0] 
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The charmonium and bottomonium ground states are thus tightly bound, with a binding 
energy AE = 2M D B — M ^> I ~ 0.2 GeV, and very small, with r « rj ~ 1 fm. What 
happens to them in a QGP? 

The effect of colour screening is that the binding becomes weaker and of shorter range. 
When the force range or the screening radius fall below the binding radius, the Q and 
Q can no longer "see" each other, and hence the bound state becomes dissociated. As 
already noted, the quarkonium dissociation points therefore determine the temperature 
and thus also the energy density of the QGP. The basic question thus is how to calculate 
the quarkonium dissociation temperatures. 

Early attempts were based on models of the heavy quark potential, essentially obtained 
from d — 1 electrodynamics. Using these in the Schrodinger equation, together with a 
crude lattice form for the screening mass, led to first predictions E] 

Tj/i> ^ Tc T x & < T c , (3) 
where T c is the deconfinement temperature. 

When lattice results for the heavy quark potential became available from finite tempera- 
ture lattice studies, these were employed in the Schrodinger equation in various forms [7 j 
- ^2] • The results eventually converged fairly well and the present status is schematically 
summarized in Table 2. 



state 


J/^{1S) 


Xc(lP) 


V(2S) 


T d /T c 


2.1 


1.2 


1.1 



Table 2: Charmonium dissociation temperatures in lattice-based potential theory 

Both previous approaches assume the validity of a two-body potential treatment at fi- 
nite temperature near a critical point. This assumption is no longer necessary if the 
quarkonium spectrum can be calculated directly in finite temperature lattice QCD. Such 
calculations have become possible in recent years and results were presented by several 
groups, first in quenched QCD [TS] - [IE] and now also in full (two-flavour) QCD [T7| ITHj. 
The present state of these results is summarized in Table 3. 



state 


J/^{IS) 


XcO-P) 


f{2S) 


T d /T c 


> 2.0 


< i.i 


? 



Table 3: Charmonium dissociation temperatures from finite temperature 

Lattice QCD calculations 

Very recently also first lattice results have been presented for bottomonium dissociation 
in quenched QCD [1211201; one finds there 

T T > 2 T c , T Xb < 1.15 T c . (4) 
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The low value reported for the Xb, which has approximately the same binding energy as 
the J/ip, remains at present quite puzzling. 

We thus find from direct finite temperature lattice studies, both in quenched and in full 
QCD, as well as in lattice-based potential work, that the J /if} and the T survive up to 
T > 2 T c , which means up to energy densities of 25 GeV/fm 3 or more. In contrast, the 
Xc and (so far only from potential studies) the if)' melt near T c , i.e., for energy densities 
in the range 0.5 - 2.0 GeV/fm 3 . It should be noted that "survival" here means that 
the corresponding signal is seen up to the temperature in question. So far, lattice QCD 
results do not yet allow a determination of the widths as function of temperature, and 
hence it is not known if even the ground states acquire a very large width with increasing 
T. Moreover, the comparison of lattice and potential theory can be carried out on a more 
detailed level than given by just the dissociation temperatures. A study of correlators in 
both approaches can thus certainly provide more insight [T2*] . 

In closing this section, we summarize the modifications in theory which have led to the 
present new theoretical view, in particular of J /if) survival in a hot QGP: 

• Earlier lattice studies had provided only the colour average of the in-medium QQ 
free energy; now it is possible to separate out the colour singlet contribution. 

• Earlier potential models had used the free energy as potential in the Schrodinger 
equation; today we can specify the colour singlet internal energy, which is more 
realistic as the relevant potential and leads to a stronger binding. 

• There now exist direct finite temperature lattice QCD studies of the in-medium 
behaviour of charmonia, allowing ab initio conclusions not based on any potential 
model, and they support a higher J/ip dissociation temperature. 

What then does this imply for quarkonium production as QGP probe in nuclear collisions? 
3. In -Medium Behaviour of Quarkonia: Phenomenology 

The modifications observed when comparing J/if> production in AA collisions to that in 
pp interactions have two distinct origins. Of primary interest is obviously the effect of 
the secondary medium produced in the collision - this is the candidate for the QGP we 
want to study. In addition, however, the presence of cold nuclear matter in target and 
projectile can also affect the production process and final rates. This ambiguity in the 
origin of any observed J /if) suppression thus has to be resolved. 

A second empirical feature to be noted is that the measured J /if) production consists of 
directly produced IS states as well as of feed-down from Xc(lP) an d V ;/ (2S) decay. In 
the previous section, we had seen that a hot QGP affects the higher excited quarkonium 
states much sooner (at lower temperatures) than the ground states. This results in another 
ambiguity in observed J /if) production - are only the higher excited states affected, or do 
all states suffer? 

An ideal solution of these problems would be to measure separately J /if), Xc and if>' pro- 
duction first in pA (or dA) collisions, to determine the effects of cold nuclear matter, 
and then measure, again separately, the production of the different states in AA colli- 
sions as function of centrality at different collision energies. While the production of the 
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ip' has been studied in both pA and AA collisions, \ c data on nuclear targets are not yet 
provided. 

Until such data become available, we resort to a more operational approach, whose basic 
features are: 

• We assume that the J/ip feed-down rates in pA and AA are the same as in pp, i.e., 
60 % direct J/ip(lS), 30 % decay of and 10 % decay of ip'(2S). 

• We specify the effects due to cold nuclear matter by a Glauber analysis of pA 
or dA experiments in terms of a % abs for i =J/ip, Xa • This a % ahs is not meant 
as a real cross-section for charmonium absorption by nucleons in the nucleus; it 
is rather used to parametrize all initial and final state nuclear effects, including 
shadowing/ antishadowing, parton energy loss as well as pre-resonance and resonance 
absorption. 

• In the analysis of AA collisions, we then use a l abs in a Glauber analysis to obtain 
the predicted form of normal J / ip suppression. This allows us to identify anomalous 
J/ip suppression as the difference between the observed production distribution 
and that expected from only normal suppression. We parametrize the anomalous 
suppression through the survival probability 

s = {dNj/dy) exp 

(dNi/dy) Glauber 

for each quarkonium state i. 

With the effects of cold nuclear matter thus accounted for, what form do we expect for 
anomalous J/ip suppression? If AA collisions indeed produce a fully equlibrated QGP, we 
should observe a sequential suppression pattern dfor J/ijj and T, with thresholds predicted 
(in terms of temperature or energy density) by finite temperature QCD 0E], |2T]-|25j. 
The resulting pattern for the J / if) is illustrated in Fig. El 



EI 



(2S)(1P) (IS) 



E(2S)E(1P) £(1S) 
Energy Density 

Figure 3: Sequential J ftp suppresssion 

Its consequences are quite clear. If, as present statistical QCD studies indicate, the direct 
J/ip{lS) survives up to about 2 T c and hence to e > 25 GeV/fm 3 , then all anomalous 
suppression observed at SPS and RHIC must be due to the dissociation of the higher 
excited states Xc and ip'. The suppression onset for these is predicted to lie around e ~ 1 
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GeV/fm 3 , and once they are gone, only the unaffected J/ipiAS) production remains. 
Hence the J/ip survival probability (under anomalous suppression) should be the same 
for central Au — Au collisions at RHIC as for central Pb — Pb collisions at the SPS. 

A further check to verify that the observed J j ip production in central collisions is indeed 
due to the unmodified survival of the directly produced IS* state is provided by its trans- 
verse momentum behaviour. Initial state parton scattering causes a broadening of the px 
distributions of charmonia j2^-[22j: the gluon from the proton projectile in pA collisions 
can scatter a number of times in the target nucleus before fusing with a target gluon to 
produce a cc. Assuming the protonic gluon to undergo a random walk through the target 
leads to 

(Pt) p a = (P 2 T ) PP + N% (6) 

for the average squared transverse momentum of the observed J/ip. Here N A specifies 
the number of collisions of the gluon before the parton fusion to cc, and So the kick it 
receives at each collision. The collision number N A can be calculated in the Glauber 
formalism; here a a b s has to be included to take into account the presence of cold nuclear 
matter, which through a reduction of J/ip production shifts the effective fusion point 
further "down-stream" [2Hj . 

In AA collisions, initial state parton scattering occurs in both target and projectile, and 
the corresponding random walk form becomes 

(p 2 t)aa = (P 2 T ) PP + N? A 5 ; (7) 

here N AA denotes the sum of the number of collisions in the target and in the projectile, 
prior to parton fusion. It can again be calculated in the Glauber scheme including a a b s - 
The crucial point now is that if the observed J/^'s in central AA collisions are due to 
undisturbed IS" production, then the centrality dependence of the px broadening is fully 
predicted by such initial state parton scattering j2Sj. In contrast, any onset of anomalous 
suppression of the J/ijj(lS) would lead to a modification of the random walk form [2%] . 

In Fig. |3] we summarize the predictions for J/ ip survival and transverse momentum be- 
haviour in AA collisions at SPS and RHIC, as they emerge from our present state of 
knowledge of statistical QCD. Included are some preliminary and some final data; for a 
discussion of the data analysis and selection, see ref. 
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Figure 4: J/ip survival and transverse momentum at SPS and RHIC 
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We conclude that the present experimental results are compatible with the present in- 
formation from statistical QCD. This was not the case previously, and such a conclusion 
can be drawn today because of several changes in our theoretical and experimental un- 
derstanding: 

• As already noted, statistical QCD presently puts the onset of direct J/ip suppression 
at energy densities beyond the RHIC range; previous onset values were much lower 
(see, e.g., ref. 

• SPS In — In data [2H] suggest an onset of anomalous suppression at e ~ 1 GeV/fm 3 ; 
previous onset values from Pb — Pb and S — U interactions were considerably higher, 
with e ~ 2 - 2.5 GeV/fm 3 (see, e.g., ref. jSD]). 

• within statistics, there is no further drop of the J/ip survival rate below 50 - 60 %, 
neither at RHIC nor at the SPS; a second drop in very central SPS Pb — Pb data 
(see, e.g., ref. jSOj) is no longer maintained. 

4. J/ip Enhancement by Regeneration 

In this section we want to consider the possibility that the medium produced in high 
energy nuclear collisions differs from the deconfmed state of matter studied in finite tem- 
perature QCD. The basic idea here is that nuclear collisions initially produce more than 
the thermally expected charm, and that this excess, if it survives, may lead to a new form 
of combinatorial charmonium production at hadronization. 

A crucial aspect in the QGP argumentation of the previous sections was that charmonia, 
once dissociated, cannot be recreated at the hadronization stage, since the abundance of 
charm quarks in an equilibrium QGP is far too low to allow this. The thermal production 
rate for a cc pair, relative to a pair of light quarks, is 

cc/qq ~ exp{-2m c /T c } ~ 3.5 x 10~ 7 , (8) 

with m c — 1.3 GeV for the charm quark mass and T c = 0.175 GeV for the transition tem- 
perature. The initial charm production in high energy hadronic interactions, however, is 
a hard non-thermal process, and the resulting rates from perturbative QCD are consid- 
erably larger. We illustrate this for pp collisions in Fig. with cc/qq = a C c/a in [3*H 132"]. 
Moreover, in AA interactions the resulting c/c production rates grow with the number 
N co u of nucleon-nucleon collisions, while the light quark production rate grows (at least 
in the present energy regime) essentially as the number N part of participant nucleons, i.e., 
much slower. At high collision energies, the initial charm abundance in AA collisions is 
thus very much higher than the thermal value. What happens to this excess in the course 
of the collision evolution? 

The basic assumption of the regeneration approach j3H| - [H5| is that the initial charm 
excess is maintained throughout the subsequent evolution, i.e., that the initial chemical 
non-equilibrium will persist up to the hadronization point. If that is the c from 

a given nucleon-nucleon collision can at hadronization combine with a c from a different 
collision ("off-diagonal" pairs) to create a J/ip. This pairing provides a new exogamous 
charmonium production mechanism, in which the c and the c in a charmonium state 
have differen parents, in contrast to the endogamous production in a pp collision. At 
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Figure 5: Thermal vs. hard charm production in pp collisions 

sufficiently high energy, this mechanism will lead to enhanced J/ip production in AA 
collisions in comparison to the scaled pp rates. When should this enhancement set in? 

In present work [33] - |3S], it is first assumed that the direct J/ip production is strongly 
suppressed for e > 3 GeV/fm 3 . This is evidently in contrast to the statistical QCD results 
discussed in the previous sections; however, we recall the caveat that the temperature de- 
pendence of the charmonium widths is so far not known. Moreover, it is of course always 
possible that the medium produced in nuclear collisions is quite different from the quark- 
gluon plasma of statistical QCD. Next, it is either assumed that the regeneration rate is 
determined by statistical combination in a QGP [33J or a specific in-medium cc recom- 
bination process is invoked, depending on the expansion geometry and the momentum 
distribution of the produced charm quarks [SI] . In general, however, if a c and a c meet 
under the right kinematic conditions, they are taken to form a J/ip. An evolution towards 
a QGP in chemical equilibrium would also require annihilation at this point. 

To account for the J/ip production rates observed at RHIC, it is assumed that the new 
exogamous production just compensates the proposed decrease of the direct endogamous 
IS production, as illustrated in Fig. El At the LHC, with much higher energy densities, 
one should then observe a J/ip enhancement relative to the rates expected from scaled pp 
results. 
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Figure 6: Regeneration of J/ip production 

We thus have to find ways of distinguishing between the two scenarios discussed here: the 
sequential suppression predicted by an equilibrium QGP or a stronger direct suppression 
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followed by a J/ip regeneration in a medium with excess charm. Fortunately the basic 
production patterns in the two cases are very different, so that one may hope for an 
eventual resolution. 

The overall J/i/j survival probability in the two cases is illustrated in Fig. [T^i. Sequential 
suppression provides a step-wise reduction: first the higher excited charmonium states 
are dissociated and thus their feed-down contribution disappears; at much higher temper- 
ature, the J/ip(lS) itself is suppressed. Both onsets are in principle predicted by lattice 
QCD calculations. In the regeneration scenario, the thermal dissociation of all "diagonal" 
J/i/j production is obtained by extrapolating SPS data to higher energy densities. The 
main prediction of the approach is therefore the increase of J / ip production with increas- 
ing energy density. Ideally, the predictions for the LHC are opposite extremes [23 | l36 t l37j. 
providing of course that here the feed-down from S-decay is properly accounted for. 




„ „ Centrality 
Energy Density 

(a) ' 0>) 

Figure 7: Sequential suppression vs. regeneration: J / ip survival (a) and p^-behaviour (b) 

The expected transverse momentum behaviour in the two cases is also quite different. In 
the region of full J/ip{lS) survival, sequential suppression predicts the normal random 
walk pattern specified through pA studies; the eventual dissociation of direct J/i/;(lS) 
states then leads to an anomalous suppression also in the average p\ [2H]- Regeneration 
alone basically removes the centrality dependence, since the different partners come from 
different collisions. It is possible to introduce some small centrality dependence [SB], but 
the random walk increase is essentially removed. The resulting behaviour is schematically 
illustrated in Fig. 03- - More generally, the quarkonium momentum distributions, whether 
transverse or longitudinal, should in the regeneration scenario be simply a convolution of 
the corresponding open charm distributions; this provides a further check |35] . 

5. Conclusions 

• In statistical QCD, the spectral analysis of quarkonia provides a well-defined way 
to determine the temperature and energy density of the QGP. 

• If nuclear collisions produce a quark-gluon plasma in equilibrium, the study of 
quarkonium production can provide a direct way to connect experiment and statis- 
tical QCD. 

• For a QGP with surviving charm excess, off-diagonal quarkonium formation by 
statistical combination may destroy this connection and instead result in enhanced 
J/tp production. 
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